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In quantum optimal control theory the success of an optimization algorithm is highly influenced by 
how the figure of merit to be optimized behaves as a function of the control field, i.e. by the control 
landscape. Constraints on the control field introduce local minima in the landscape -false traps- 
which might prevent an efficient solution of the optimal control problem. Rabitz et al. [Science 
303, 1998 (2004)] showed that local minima occur only rarely for unconstrained optimization. Here, 
we extend this result to the case of bandwidth-limited control pulses showing that in this case one 
can eliminate the false traps arising from the constraint. Based on this theoretical understanding, 
we modify the Chopped Random Basis (CRAB) optimal control algorithm and show that this 
development exploits the advantages of both (unconstrained) gradient algorithms and of truncated 
basis methods, allowing to always follow the gradient of the unconstrained landscape by bandwidth- 
limited control functions. We study the effects of additional constraints and show that for reasonable 
constraints the convergence properties are still maintained. Finally, we numerically show that 
this approach saturates the theoretical bound on the minimal bandwidth of the control needed 
to optimally drive the system. 

PACS numbers: 02.30.Yy,02.60.Pn,03.65.Aa,03.67.-a 


The ability of achieving a desired transformation of a 
quantum system lies at the heart of the success of ex¬ 
periments in cold atoms mm, quantum optics ma, 
condensed matter [6], and in quantum technologies [7j. 
Together with the fast development of these fields and 
the increasing complexity of the experiments, to develop 
efficient protocols it is often necessary to automatize the 
optimization process if not the whole development of the 
experimental sequence [5]. One possible way to perform 
such an optimization is by means of quantum optimal 
control, an approach that has proven to be very success¬ 
ful in solving this class of problems [5]. 

The success of gradient methods to find global solu¬ 
tions in quantum control problems [TUHT^ is largely due 
to the fact that local minima are very rare in the control 
landscapes of a large class of systems UMI]. However, 
to face the new challenges posed by the recent advance¬ 
ments in quantum science, gradient methods might be 
not the best option since to numerically calculate the 
gradient of the control objective might be quite inef¬ 
ficient. Moreover in an increasing number of interest¬ 
ing applications, the control objective does not allow for 
an analytical calculation of the gradient [22l [23]. The 
CRAB optimal control algorithm operates with an ex¬ 
pansion of the control field onto a truncated basis and 
a direct optimization of the coefficients of the expansion 
by means of gradient-free minimization [241 ES] ■ These 
characteristics allow for the solution of optimal control 
problems involving many-body quantum systems |2^l26j . 
as well as in the presence of long-range interactions 
m, bandwidth-limited control according to experimen¬ 
tal constraints pHlldOj . and highly nonlinear functionals 


where the gradient can not be calculated [53]. How¬ 
ever, the CRAB optimization does not necessarily fulfill 
the condition under which local minima in the optimiza¬ 
tion occur only rarely as it is by construction bandwidth- 

limited [IMS]. 

In this paper we extend the results presented by Rabitz 
and coworkers m and show that bandwidth-limited op¬ 
timal control can be made virtually free of local minima 
as in the case of unconstrained control. Moreover, we nu¬ 
merically show that the global minima one always reaches 
correspond to an optimal solution if the bandwidth sat¬ 
isfy the theoretical bound given in Reference m- In par¬ 
ticular, we present an extension of the CRAB algorithm - 
the “dressed CRAB” (dCRAB) - that keeps the benefits 
of the original algorithm and comes with the additional 
property of guaranteed convergence to the global opti¬ 
mum in the cases where this is guaranteed also for gradi¬ 
ent methods. We test the two versions of the algorithm 
by optimizing a state transfer for different instances of 
a random spin Hamiltonian. Being M the dimension 
of the Hilbert space the problem is defined on, for the 
standard CRAB algorithm false traps occur when the 
heuristic “2M—1-rule” for the number of required control 
coefficients is violated |33|. On the contrary, using the 
dCRAB presented hereafter, all false traps are removed 
regardless of the number of control coefficients, resulting 
also in a faster convergence to the global optimum. We 
present a theoretical explanation supporting these find¬ 
ings and showing that one can construct a set of random 
basis functions which follow the instantaneous gradient 
of the control landscape. Finally, we examine the behav¬ 
ior of dCRAB in the presence of additional constraints on 
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the bandwidth and amplitude of the pulse, and show that 
the algorithm is well-behaving as false traps appear only 
in the presence of strong constraints. In the following, 
for simplicity we will focus on state to state transfer of 
pure states, however the theoretical arguments are valid 
in the general scenario of unitary gate generation and 
mixed-states optimal control. 

The structure of this paper is as follows: we first review 
the standard CRAB algorithm and introduce its exten¬ 
sion, dCRAB [33]. In section II we present the theoret¬ 
ical basis of dCRAB, while in section III we show that 
dCRAB follows the instantaneous gradient of the land¬ 
scape. Finally, in Section IV we apply both CRAB and 
dCRAB to control problems with random spin Hamilto¬ 
nians and evaluate their performance with respect to suc¬ 
cess probability, computational effort and behavior under 
constraints. 


I. FROM CRAB TO DRESSED CRAB 

In the following, we consider the optimization of a 
state-to-state transfer problem where the figure of merit 
is the overlap of the final state after the evolution with 
the target state 

FmT))) = \{(:mT))\\ (i) 

Here |C) is the target state and the time evolution is given 
by the Schrodinger equation 

\m) = \o ( 2 ) 

with initial state |^) and control function f{t) is deter¬ 
mined by the optimization algorithm in order to maxi¬ 
mize the figure of merit. 

The CRAB algorithm builds on the fact that in 

most scenarios the resources available to solve an optimal 
control problem - such as time, energy and bandwidth 
- are limited: in particular, as the set of practically ac¬ 
cessible wave functions is usually limited, one can show 
that also the control bandwidth of the control field can 
be upper bounded |3T]. This bound has very important 
consequences as optimal control problems can be prac¬ 
tically solved by exploring a small subset of the a priori 
infinite dimensional search space of functions. That is, 
one can expand the control field in a truncated basis 


mim 




FIG. 1. (Color online) Schematic view on the control land¬ 
scape J(/) = F{\ip{T))). While F has a strict global max¬ 
imum point at 1^) (left panel), this state can be reached by 
different control functions fi and /2 corresponding to multiple 
global maximum points in the landscape J (right panel). 


standard choice for the basis functions fi(t) are trigono¬ 
metric functions, often multiplied by a shape function 
l/ 7 (t) that fixes the pulse boundary conditions and a 
guess pulse [331135] . 

Due to the restriction of the search basis to Nc di¬ 
mensions given by the CRAB expansion in Eq. (|^, the 
algorithm might converge to a non-optimal fixed point, 
i.e. the algorithm is trapped in a local minimum arising 
due to the constraint - a so-called false trap. To over¬ 
come this problem, and escape from these false traps, 
we show in the next section that one can start from the 
non-optimal fixed point a new CRAB optimization with 
a new random basis and new coefficients. This is done in 
an iterative way so that in the j-th super-iteration one 
optimizes the coefficients of 


Nc 

+ (4) 

i=l 

where {t) are new randomly chosen basis functions, for 
example sine or cosine functions with random frequencies 
within some interval [0, Wmax]- As a consequence, in each 
super-iteration the old pulse is dressed with new search 
directions and we call this procedure dressed Chopped 
Random Basis (dCRAB) algorithm [33]. This updating 
of the search directions can also be understood as an ex¬ 
tension of Powell’s method [84j to an infinite dimensional 
search space. In the following sections we give a theoret¬ 
ical explanation why this is a substantial improvement 
of the algorithm by analyzing how it influences the con¬ 
strained control landscape and we demonstrate it numer¬ 
ically by applying it to a model with typical properties. 


II. CONTROL LANDSCAPES AND CRAB 


Nc 

fit) = (3) 

z=l 

and the optimization can then be performed on this sub¬ 
space of small dimension, resulting in an optimal set of 
coefficients Ci,i = 1, • • • ,Nc- The optimization can be 
performed by standard tools, e.g. by the Nelder-Mead 
simplex algorithm that does not rely on gradients (3^ . A 


In this and in the following section, we review the the¬ 
ory of control landscapes and how it can enlighten the 
reasons why optimal control algorithms converge to the 
optimal solution or become trapped [31 [T3|. We first 
review a general perturbative analyses of control land¬ 
scapes presented in specifically its gradients and crit¬ 
ical points, valid independently from the particular opti¬ 
mization algorithm employed to find the optimal driving 
field; and then we specify it to the CRAB approach we 
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are focusing here. In particular, we focus on why CRAB 
can be trapped whereas dCRAB - like gradient methods 
- cannot be trapped, as we will show also in the next 
section by means of some numerical examples. 

Given a control problem with control f(t), the control 
landscape is the functional J(/) with 

J(/) = F(|V^(T))) (5) 

where \'ijj{T)) is the final state resulting from time evo¬ 
lution of the initial state |^) with the given control f{t) 
and is a figure of merit quantifying the quality of the 
process mm- In other words, the control landscape is 
the functional expressing the figure of merit of the pro¬ 
cess as a function of the control field. In the following we 
use the state fidelity given by Eq. Q as a figure of merit. 
Notice that the control landscape can be seen either as 
a function of the final state \tjj{T)) or as a functional of 
the control field f{t): in the first case the function has by 
definition a single optimal point {\ip(T)) = |^), up to a 
global phase), while in the latter there might exist differ¬ 
ent control pulses leading to the same maximum value, 
as sketched in Fig. 

Optimization usually leads to so-called critical points 
of the landscape, that is to pulse shapes fulfilling the 
condition 


5J=(VF(V^(r))|(5V'(T)) = 0 ysf, (6) 

i.e. a vanishing variation of the functional J for a varia¬ 
tion of the control /. Using the chain rule one can show 
that this variation consists of two parts: the gradient of 
the fidelity as a function of the final state WF/5ij}{T), 
and the variation of the final state as a result of the vari¬ 
ation of the control S'ip{T)/6f. For all common choices 
of the figure of merit a vanishing gradient WF{tjj{T)) = 0 
corresponds to the global maximum, global minimum or 
a saddle point laiHin]; more specifically if F{\^|;{T))) 
is the state overlap (eq. 0 ) there are no saddle points 
since VF = 0 corresponds to F = 0 (the global mini¬ 
mum) or F = 1 (the global maximum). Less intuitive 
is the role played by the second term, which we analyze 
in the following, to understand under which conditions 
SJ = 0 implies VF(i/;(F)) = 0. 

Critical points can be classified as singular or regu¬ 
lar critical points, where regular means that for every 
|5'0(T)) in the Hilbert space of the problem there is a 
change in the control Sf that generates it, while on the 
contrary for singular points this is not the case [smiin]. 
As a direct consequence, for regular critical points we 
have 


SJ = {\/F{jp{T))\S't{j{T)) = 0 V (5/ 
^SJ= {WF{'ilj{T))\d^{T)) = 0 V S-tpiT) 
^ VF(V'(T)) = 0 


(7) 


point is not a trap. That is, all traps have to be sin¬ 
gular critical points (although the reverse is not neces¬ 
sarily true). However, at least for controllable systems, 
all known traps occur at constant control fields nmis], 
often at / = 0, and numerical evidence hints to the fact 
that the singular points at non-constant control seem to 
be no traps [TMT5] . That is, it is commonly accepted 
that unconstrained control landscapes for all practical 
purposes have no traps, although a rigorous proof has 
been given only for two specific situations, the Landau- 
Zener system and the transmission of a wave package 
over a potential barrier m- 

Now we discuss the implications of these results on 
the CRAB and dCRAB optimal control success rate. As 
a starting point, we recall that the previous statements 
assume that the control / is a general function in L 2 
whereas for CRAB we are much more restricted by the 
expansion into the truncated basis. This means that one 
might find false traps in the CRAB control landscape, i.e. 
encounter points where SJ vanishes for all variations Sf 
allowed by the truncated basis expansion of Eq. (^, but 
not for all variations of the unconstrained control space: 

No 

(5J = 0 V <5/ = ^ fi{t)Sci but 

i=l 

^^f'f + Sf€L2,SJf=0. ( 8 ) 

These points are called false traps [32] as they arise only 
artificially from the choice of the basis and their influence 
can hinder convergence of the algorithm (32] [3^ . 

Here, we show how these false traps for CRAB are 
removed by the super-iterations of dCRAB: From equa¬ 
tions and ^ together with 

(VF(|^(T)))|.) = 2Re((V'(F)|C)(C|-)) (9) 

it follows that, for a false trap the gradient has to be 
(VF('(/)(T))| •) = Re(0T|-) for some non-zero vector 
\(j)T)- A perturbative treatment yields [17]: 

\Si;{T)) = -iU{T) r uHt)H,um)Sf{t)dt . (10) 

Jo 

This results in the expression for the overlap of the state 
update and the gradient 

SJ = Re{(j)T\Stlj{T))= [ k{t)Sf{t)dt (11) 

Jo 

k(t) = -Im {(PT\U{T)U^{t)H^U{t)\Ci (12) 

where fc is a continuous function and k d since 
fc = 0 would violate Eq. If we now choose 6f{t) = 
sm{uJrt)Sc with a new random frequency ujr we get 

k{t)Sf{t)dt ^ 0 (13) 



as in this case the gradient of the fidelity is orthogonal almost surely (i.e. the integral vanishes only on a null 
to the whole Hilbert space. This means that a regular set of the probability measure). If we now perturb our 
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control field / by this new frequency contribution (this 
is the new super-iteration of dCRAB) we find 

5J={VF{i;{T))\5i,(T))^Q (14) 

and thus we have removed the false trap. In conclusion, 
the recipe to escape from a false trap is simply to add 
a new random frequency term to the CRAB expansion 
once we are in the false trap. This can be done also 
without increasing the total number of coefficients since 
when one is at the bottom of a false trap for a given set 
of basis functions and coefficients there is no more use in 
varying the old coefficients, and they can be kept at their 
value. We thus use the pulse / leading to the false trap 
as a guess pulse for the next super-iteration of dCRAB 
(as explained in the previous section, see Eq. @) that 
opens the false trap. 

III. INSTANTANEOUS BASIS FUNCTIONS 
FOR PULSE UPDATE 


out of the whole tangential space of the unity sphere of 
(i.e. regularity of the control landscape at the point 
of consideration which is the general case for trap free 
landscapes). The state update generated by each pulse 
update in a perturbative regime reads m- 

\d^niT)) = -iU{T) [ uHt)HiU{t)\OSfn{t)dt. (19) 
Jo 

We now show that for every vector \S(j>) out of this 
unit sphere there is a set of real coefficients a„ so that 
ctnlStpniT)). In Order to do so we identify 
\6(j)) and |<5i/'rt(T)) as elements of a 2M — 1 dimensional 
real vector space with scalar product Re (ulrc) and show 
that the \Sipn{T)) are a basis of this space. Let us thus 
consider = span{\Stf}i{T)),... ,\Stf}n{T))}. We will 
show by induction that has dimension n. For n = 1 
this is trivial. If now has dimension n we have to 
show that \SilJn+i{T)) is not orthogonal to P^. Let thus 
be u e Pn- The overlap of |5'0„+i(T)) and |u) is 


In this section we show that by choosing enough new 
frequencies in a single super-iteration of dCRAB we can 
follow approximately the instantaneous gradient of the 
control landscape with a bandwidth-limited pulse up¬ 
date. In the previous section we chose a new basis func¬ 
tion so that its scalar product with 

k{t) = -Im {(l)T\U{T)U\t)HiUitm (15) 


would be finite. It is worth noting that this function 
k(t) is exactly the update direction used in gradient al¬ 
gorithms with the time-evolved initial state 

\m) = u{tm ( 16 ) 


and the adjoint state 

\x{t)) = U{t)UHT)\VF{^{T))). 


(17) 


Together with Hi = ^ ^ ^g^. 


known 


k{t) = -Im {x{t) 


Sf 


5H, , 




(18) 


To follow the gradient of the control landscape only a 
limited freedom for the pulse update is needed. Hsieh et 
al. [37] derived for a state to state transfer in an M-level 
system a set of 2M — 2 functions (for each iteration step) 
that span the pulse update k{t). This simply reflects the 
fact that the state update generated by the pulse update 
6f, i.e. \5'tp{T)), is a M dimensional complex vector (2M 
real coefficients) with norm I and irrelevant global phase 
(thus 2M — 2 real coefficients). 

We can generate something very similar by choosing 
2M — I random frequencies (we do not fix the global 
phase here) and corresponding pulse updates 5fn{t) = 
sin(a;„t)(5c„ (n = I,... 2M — I) with real coefficients c„. 
Let us here assume that L -2 pulses can generate \5'ij){T)) 


Re{v\6tpn+i(T)) = / l{t)Sfn+i{t)dt (20) 

Jo 

l{t) = -lm{v\U{T)UHt)HiUm) 

The kernel l{t) is continuous and ii I ^ 0 the overlap is 
nonzero almost surely. However, Z = 0 |u) = a\'ijj{T)) 

{a G M) [3Sj and thus ju) is not in the tangential space. 
This proves that the jSi/^niP)) as span the whole 2M — 
1 dimensional real vector space, or in other words the 
|^■^/'n(T)) with real coefficients span the whole tangential 
space of the unit sphere of (almost surely). 

We can also orthogonalize the state updates in a Gram 
Schmidt way to obtain the orthogonal basis \5\pn{T)). 
This translates into a change of the pulse updates via 

n—1 

5Ut) = suit) - Y,{SMT)\6MT))Sp{t ), (21) 

k^l 

with \S'ijjn{T)) generated by Sfn{t). Each direction \5(j)) 
is then a linear superposition of the basis states 

2M-1 

\6cl>)= Y, c^nlSMT)) (22) 

n—1 

and can thus be generated by a superposition of the cor¬ 
responding pulses 

2M—1 

\S^)=-1 / U{T)U^{t)HiU{t)\0 Y C^nSfn{t)dt. 

"'o 

(23) 

In particular this proves that with 2M — I random in¬ 
stantaneous basis functions {Sfn{t) oder Sfn{t) for the 
“orthogonal” basis) we can follow the instantaneous gra¬ 
dient of the control landscape by suitable coefficients. 
We stress that, unlike in the case of a standard gradient 
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method, this is done by bandwidth-limited pulses which 
can be comfortably adapted to typical experimental con¬ 
straints such as bandwidth-limited control electronics. 

This result demonstrates that after a change of the 
truncated basis in the function space the search in the 
state space is quasi-local. We can use this in the nu¬ 
merical part by initializing a small simplex after the ba¬ 
sis change to ensure local search. The orthogonalization 
procedure could also help to find a better basis for the 
CRAB (re)start at the cost of M{M — l)/2 additional 
function evaluations. 


IV. NUMERICAL EXPERIMENT 

In this section we show that the dCRAB algorithm is 
indeed capable of escaping from false traps in random 
optimal control problems of increasing complexity. We 
consider a spin Hamiltonian 

N N-1 

H = J2 + Aaf + fit) '"f+i (24) 

1=1 i=l 

with random coefficients Oi, /3i S [ 0 , 1 ] and a control field 
f{t) tuning the interaction. The random coefficients lift 
the symmetry of the system and make it controllable [33] ■ 

We investigate the control resources needed to drive 
the system from a random initial state |^) to a ran¬ 
dom target state |C) within a fixed time interval [0,T]. 
We measure the fidelity of this state-to-state transfer by 
Eq. Q. We investigate different system sizes and differ¬ 
ent operation time intervals, analyze the occurrence of 
false traps in the landscape by several trials of optimiza¬ 
tion of the fidelity F for different random instances of 
1 ^) and 1 C) (uniformly distributed over the unit sphere), 
and measure the success of an optimization by a certain 
threshold ry = 10“^ for the residual error e = 1 — F. 
We count an optimization trial as a success if £ < 77 
after convergence of the optimization or after a certain 
maximum number of function evaluations. The success 
rate of the optimization depends on the presence or ab¬ 
sence of false traps in the landscape. To detect them 
we optimize ten different sets of random instances of the 
Hamiltionian and of the initial and goal state; for each 
instance we use ten different random starting points and 
random sets of frequencies to perform the CRAB opti¬ 
mization. The frequencies of the basis functions (both 
for CRAB and dCRAB) are chosen randomly out of an 
interval [0,u}max]- we choose {2tt) = 8,20,40 for 

N = 2, 3,4 respectively. 

Fig. [^reports the success probability pc for the stan¬ 
dard CRAB optimization as a function of the number 
Nc of coefficients of the truncated basis expansion in 
Eq. (§: as can be clearly seen, for Nc large enough no 
false traps are present, resulting in a success probability 
of one. Here, by large enough we mean the empirical 
‘Wc = 2 • 2 ^ — 2 ” rule for unconstrained optimization, 
where the number of real coefficients equals the number 



FIG. 2. (Color online) Success probability pc as a function 
of the number Nc of basis functions in the CRAB expansion. 
As we increase the size of the function space false traps are 
removed. The symbols are the mean values over 10 different 
random pairs of Hnal and initial states |^) and |C) with 10 
different starting points and frequencies each. The error bars 
show the standard deviation over the different final and initial 
states. The total time was set to T = Ott, IOtt, IOtt and the 
allowed bandwidth to ra;niax/(27r) = 8,20,40 for N = 2,3,4 
qubits respectively. An optimization was counted as success 
when the residual error e was smaller than 77 = 10~^. The 
black lines indicate the empirical “2 • 2^ — 2 ” rule for the 
required number of coefficients Nc- 


of independent real entries in the state vector [33]. The 
correspondent analysis for the dCRAB approach is re¬ 
ported in Fig.[^ the success probability pd is always one 
regardless of Nq, false traps are avoided and the fidelity 
always exceeds the threshold. 

Despite the striking difference in terms of success prob¬ 
ability, an important benchmark for any optimization 
method is given by the computational effort required to 
arrive at the optimal solution. Here we focus on the num¬ 
ber of function evaluation needed to achieve the global 
optimum as this is practically the only difference between 
the two methods: Fig. [^ shows the number of function 
evaluations Uf required by the two methods to exceed 
the threshold F = 1 — ry as a function of the coefficients 
Nc in the CRAB expansion (in the case of dCRAB the 
coefficients of a single super-iteration). All points consist 
of the average number of function evaluations of the suc¬ 
cessful runs divided by the respective success probability, 
that is the number of function evaluations that on aver¬ 
age one has to do to solve the optimal control problem 
using one of the two methods. Note that the minimal 
effort does not follow the “2 • 2^ — 2” rule of Fig. for 
guaranteed convergence. For CRAB the computational 
effort heavily depends on Nc, which can be problematic 
as the best choice of Nc is not known in advance. For 
dCRAB instead there is only a minor dependance on Nc 
in the order of magnitude of the error bars. Furthermore, 
Fig.j^shows that even with the best choice of Nc CRAB 
can not beat the performance of dCRAB. 
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FIG. 3. (Color online) Success probability pd as a function of 
the number Nc of basis functions in a single call of CRAB 
within the dCRAB super-iterations. Independently of Nc 
no false traps are encountered and pd = 1. The total time, 
allowed bandwidth and error threshold are as in Fig. 



FIG. 4. (Color online) Number of function evaluations n/ 
for dCRAB (red) and original CRAB (black) as a function of 
the number Nc of basis functions involved in a single call of 
the respective algorithm. The error bars show the logarithmic 
standard deviation. Optimization was stopped when the error 
£ crossed the threshold p — 10~^. The total time was T = 
IGtt, the allowed bandwidth TcOmaxI (Stt) = 40 and the system 
size A = 4 qubits. Similar results are obtained for other 
choices of the parameters. 


A. Constrained optimization 


Finally, we study how the convergence properties 
change in the presence of additional constraints - like 
limited fluence and bandwidth typically present in ex¬ 
perimental setups - that violate the hypothesis of the 
analysis of control landscapes presented up to now. In¬ 
deed, in this scenario, false traps might be present which 
might change performance and convergence speed of the 
algorithms. We consider separately two different kinds of 
constraints: bandwidth-limited control and limited pulse 
height. 

In the first case, we observe that Eq. (13) still holds 
even if the new random frequency ojr is chosen only 
within the limited bandwidth interval [0,a;max]- We can 
then study the performance of the optimization as a func¬ 
tion of Wmax, as done in the unconstrained case and com- 
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FIG. 5. (Color online) Infidelity £ = 1 — F as a func¬ 
tion of the limited bandwidth iJmax for |^) = |0000) and 
1^) = in the system of A = 4 qubits and total time 

T = IGtt with fixed random coefficients Oi, in the Hamilto¬ 
nian (Eq. ( |24[ )). The grey diamonds show the infidelity £ for 
CRAB (Ac = 40), while the green circles show it for dCRAB 
(Ac = max{2a;maxr/27r,40}). The black dashed line indi¬ 
cates the error threshold of = 10“^. The error bars report 
the logarithmic standard deviation. 


pare the CRAB and dCRAB approaches. The results 
are reported in Fig. where we show the optimal in¬ 
fidelity 1 — F reached from ten independent runs as a 
function of Umax- The optimal control problem is to per¬ 
form the state transfer from |^) = |0000) to |C) = 
given the Hamiltonian of Eq. (24) for CRAB (grey di¬ 
amonds) and dCRAB (green dots). One can see that 
for ujmax > 16 • 2tt/T dCRAB succeeds with probability 
one, as all instances reached the optimal fidelity. Notice 
that this is less than half the bandwidth of the previ¬ 
ous results reported in Eigs. and In addition, in 
a small intermediate regime around uJmax = 14 • 21 : jT^ 
some optimizations succeed while others fail indicating 
the presence of false traps (note that the graph shows 
just the mean value of the infidelity e and the standard 
deviation); while for uJmax < 10 • 27r/T the final state 
cannot be reached anymore, indicating that the band¬ 
width is too small to achieve the desired result. In the 
case of CRAB, the three regimes are shifted toward larger 
frequencies. The lower bound observed in both cases is 
in agreement with an information theoretical argument 
given in m- to achieve full control over the system the 
inequality UmaxT > D has to be fulfilled, where D = 32 
is the dimension of the state space. This inequality ba¬ 
sically says that the control has to contain enough in¬ 
formation to distinguish the target state from the other 
reachable states and it yields uJmax ^ 5 • 27r/T, the value 
of the bandwidth where the infidelity in Eig. [^starts to 
drop indicating that control over the system starts to be 
effective. 


We then perform a similar analysis for pulse height 
limited control with dCRAB, where we study two sce¬ 
narios to include such a constraint. We first introduce a 
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FIG. 6. (Color online) Infidelity e = 1 — F as a function of 
the limited pulse height fmax for the transfer from |^) = |00) 
to |(() = |11) in the 2-qubit system with constrained dCRAB. 
The black line indicates the error threshold of ?7 = 10~^. The 
orange triangles are obtained with Lagrange multipliers (see 
Eq. ( |25[ )), while the brown circles are obtained with a cut-off 
at fmax (see Eq. The error bars indicate the standard 

deviation over 10 different starting points of the optimization. 
The orange triangles have also errorbars in /max since the 
Lagrange multiplier is not a hard wall and instead we plot 
the maximum absolute value the optimal pulse takes in each 
realization. 

smooth constraint as usually done, that is a penalty on 
the pulse height so that the control objective becomes 

J = F — Amax |/(t)| . (25) 

As a second alternative, we limit the pulse height by a 
hard wall constraint, by using the update formula 

Nc 

= f ( 26 ) 

fUt) = I ^ (27) 

'(sgn(/'^(t))/i„ax otherwise. 

The results of these two procedures are reported in Fig.[^ 
where we plot the inhdelity 1—F as a function of the max¬ 
imal pulse height fmax- Clearly the optimization works 


better with the hard boundaries than with the Lagrange 
multiplier, as the fidelity threshold 1 — ?7 can be exceeded 
for about three times weaker pulses. This difference can 
be understood by the fact that the hard wall introduces 
pulses of higher bandwidth. Compared to the uncon¬ 
strained system we can decrease the maximal value of 
the pulse by a factor of fifteen, while keeping 100 % suc¬ 
cess probability. For even smaller cut-off /max the small 
error bars indicate that optimization failure is most prob¬ 
ably more due to a loss of controlability than due to false 
traps. 

V. CONCLUSIONS 

We have generalized the results presented in m to the 
case of bandwidth-limited control pulses, under the con¬ 
dition that the bandwidth satisfy the theoretical bound 
introduced in m- 

Thanks to this theoretical result, we have modified 
the CRAB optimization algorithm to efficiently combine 
the advantages of gradient methods with those of trun¬ 
cated basis methods. We showed that it is possible to 
exploit both the guaranteed convergence to the global 
optimum that gradient algorithms exhibit in the frequent 
case in which the kinematic (F(|^(F)))) and the dynam¬ 
ical (J(/)) landscapes are equivalent and the numerical 
gradient-free truncated basis approach. Moreover, we 
showed for two typical constraints, namely the limited 
fluence and the limited bandwidth constraints below the 
theoretical bound, that some of the convergence proper¬ 
ties survive, if the constraints are carefully implemented 
in the optimization procedure. 

We expect that the presented results will allow to 
tackle in the near future both theoretically and exper¬ 
imentally even more complex many-body problems as 
done so far, as well as a broader variety of control ob¬ 
jectives and constraints. 
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